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Abstract 

This is a summary of a central argument in recent review articles 
by the author (Physica A 313, 252 (2002), Annals of Physics 303, 226 
(2003), and Rev. Mod. Phys, July 2003). An effective field theory is 
derived for the low energy spin singlet excitations in a paramagnetic 
Mott insulator with collinear spin correlations. 



1 INTRODUCTION 

In a recent article (intended for an audience of experimentalists), the au- 
thor has reviewed arguments that many aspects of the physics of the cuprate 
superconductors can be understood by using their proximity to paramagnetic 
Mott insulators. Further, a distinction was made between Mott insulators 
with collinear and non-collinear spin correlations, and it was argued that cur- 
rent experimental evidence suggests that we need only consider the collinear 
class. A phenomenological description of the ground states and excitations 
of these classes of Mott insulators was provided, along with a discussion of 
their experimental implications. A more technical discussion (intended for 
theorists) of such insulators, along with a description of the effective field 
theories which describe their low energy properties appears in Ref. [21 E]- 
Here, we briefly recall the derivation and properties of the effective field the- 
ory of Mott insulators with collinear spin correlations, which is expressed in 
terms of a compact U(l) gauge field. The non-collinear class leads naturally 
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to a Z2 gauge theory, but we will not consider it here. The reader is referred 
to these previous reviews [H 121 El for complete citations to the hterature. 



2 Compact U(l) gauge theory of Mott insu- 
lators 

We focus on Mott insulators on a d dimensional bipartite lattice of sites j. 
The SU(2) spin operator on site j at imaginary time r can be written as 

Sj-(r) =r/j.Sn(rj,r); (1) 

here rij = ±1 on the two sublattices, Vj is the spatial co-ordinate of site j, 
n is a unit length vector in spin space, and S is the (integer or half-odd- 
integer) angular momentum of each spin. The antiferromagnetic exchange 
interaction between near neighbor spins implies that n(rj, r) will be a slowly 
varying function of its spacetime arguments. A standard analysis of the 
coherent state path integral of over the SU(2) Sj spins shows that the low 
energy quantum fluctuations are described by the following partition function 

Z = / V.(r, r) - 1) exp [-,S ^ * / A(n(.„ r)) 



-!- [ d'rdr ({ariif + c''(V,n)2) 

'^QC J 



(2) 



where c is the spin-wave velocity, and g is a coupling constant which controls 
the strength of the quantum fluctuations. Excluding the first Berry phase 
term, this is the action of the 0(3) non-linear sigma model in ci-|- 1 spacetime 
dimensions. Here we are primarily interested in the consequences of the Berry 
phases: ^T-(n(r))(ir is defined to be the oriented area of the spherical triangle 
defined by n(r), n(r + dr), and an arbitrary reference point Uq (which is 
usually chosen to be the north pole). 

The theory Q can be considered to be the "minimal model" of antiferro- 
magnets. In dimensions 0? > 1 it has at least two phases: at small g there is 
the conventional magnetically ordered "Neel" phase with (n) 7^ 0, while at 
large g there is a "quantum disordered" paramagnetic phase which preserves 
spin rotation invariance with (n) = 0. We are especially interested here in 
the nature of this paramagnetic state. In this section, we will manipulate 
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Z in this large g regime, and derive an alternative formulation which allows 
easier computation of the integral over the Berry phases. 

The key to an analysis of the large g regime is a better understanding 
of the nature of Ar- We will see that At behaves in many respects like the 
time-component of a compact U(l) gauge field, and indeed, this accounts for 
the suggestive notation. All physical results should be independent of the 
choice of the reference point no, and it is easy to see by drawing triangles on 
the surface of a sphere that changes in no amount to gauge transformations 
of At- If we change no to n'g, then the resulting A!^ is related to At by 

= A - 9.0(r) (3) 

where 0(r) measures the oriented area of the spherical triangle defined by 
n(r), no, and nj,. Furthermore, as we will discuss more completely below, 
the area of any spherical triangle is uncertain modulo 47r, and this accounts 
for the 'compactness' of the U(l) gauge theory. 

We proceed with our analysis of Z. First, we discretize the gradient terms 
of the 0(3) sigma model. We will limit our considerations here to antifer- 
romagnets on d dimensional cubic lattices, but similar considerations apply 
to other bipartite lattices. We also discretize the imaginary time direction, 
and (by a slight abuse of notation) use the same index j to refer to the sites 
of a (i + 1 dimensional cubic lattice in spacetime. On such a lattice we can 
rewrite Q as 

Z = dn^S (n^ - 1) exp ^ n^^ ■ Uj+f, -iSj^ Vj^jr^ ' (4) 

where the sum over fi extends over the d + 1 spacetime directions. We have 
also dropped unimportant factors of the lattice spacing and the spin-wave 
velocity in (jD). 

As noted above, we are especially interested here in the large g regime 
where there are strong fluctuations of the n^. There are strong cancellations 
from the Berry phases between different spin configurations in this regime, 
and so the second term in Z has to be treated with great care. We will 
do this by promoting the field Aj^ to an independent degree of freedom, 
while integrating out the n^. Notice that we have now introduced all d -|- 1 
components of the compact U(l) gauge field with the index fi, while only 
the = T component appears explicitly in (0)). The remaining components 
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appear naturally as suitable degrees of freedom when we integrate the 
out. Formally, the integration over the can be done by introducing new 
'dummy' variables Aj^ and rewriting Q by introducing factors of unity on 
each link; this leads to 

Z = In dA,, exp I -t2S J2 VjA.r J f U dn,5{n^ - 1)5(^-^/2 - A,,) 



X exp 



j n^^J>exp I -Sa^Aj^) - i2S^r]jAjr 



(5) 



In the first expression, if the integral over the Aj^ is performed first, we triv- 
ially return to (jU; however, in the second expression we perform the integral 
over the rij variables first, at the cost of introducing an unknown effective 
action Sa for the Aj^. In principle, evaluation of Sa may be performed order- 
by-order in a "high temperature" expansion in Xj g: we match correlators of 
the Aj^ flux with those of the Ajp, flux evaluated in the integral over the 
rij with positive weights determined only by the Xjg term in Q. Rather 
than undertaking this laborious calculation, we can guess essential features 
of the effective action Sa from some general constraints. First, correlations 
in the decay exponentially rapidly for large g (with a correlation length 
~ l/ln((7)), and so Sa should be local. Second, it should be invariant under 
the lattice form of the gauge transformation (jS)) 

4^ = A^^ - A^0,/2 (6) 

associated with the change in the reference point on the unit sphere from Hq 
to iIq, with (^j equal to the area of the spherical triangle formed by n^, no 
and IIq. Finally the area of any triangle on the sphere is uncertain modulo 
47r and so the effective action should be invariant under 

Aj^ ^ Aj^ + 27r. (7) 

The simplest local action which is invariant under © and (|7j) is that of 
compact U(l) quantum electrodynamics and so we have 



(8) 
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for large g; comparison with the large g expansion shows that the coupling 
~ g"^. In (jSI), is the discrete lattice derivative along the fi direction, 

and the sum over □ extends over all plaquettes of the d + 1 dimensional cubic 

lattice — both notations are standard in the lattice gauge theory hterature. 
The first term in the action (jHJ is, of course, the standard 'Maxwell' term 

of a compact U(l) gauge field. In this language, the Berry phase has the 

interpretation of a J J^A^ coupling to a fixed matter field with 'current' 
= 2SSf^r- This corresponds to static matter with charges ±25" on the two 

sublattices. It is this matter field which will crucially control the nature of 

the ground state. 

The remaining analysis of Z depends upon the spatial dimensionality 
d. In (i = 1, a dual model of (jSj) is solvable, and the results are in complete 
accord with those obtained earlier by Bethe ansatz and bosonization analyses 
of spin chains. We will consider the d = 2 case in the section below. There 
has been relatively little discussion of the d = 3 case (which exhibits both 
confining and deconfining phases of the gauge theory), and this remains an 
important avenue for future research. 

3 Duality mapping in d = 2 



As is standard in duality mappings, we first rewrite the partition function in 
2 + 1 spacetime dimensions by replacing the cosine interaction in (jH)) by a 
Villain sum over periodic Gaussians: 




(9) 

where e^j^A is the total antisymmetric tensor in three dimensions, and the g^-^ 
are integers on the links of the dual cubic lattice, which pierce the plaquettes 
of the direct lattice. Throughout this subsection we will use the index J to 
refer to sites of this dual lattice, while j refers to the direct lattice on sites 
on which the spins are located. 

We will now perform a series of exact manipulations on Q which will 
lead to a dual interface model [U E]. This dual model has only positive 
weights — this fact, of course, makes it much more amenable to a standard 
statistical analysis. This first step in the duality transformation is to rewrite 
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X i'^ X 
i+* X * i+* X 



Figure 1: Specification of tlie non-zero values of the fixed field a?^. The 
circles are the sites of the direct lattice, j, while the crosses are the sites of 
the dual lattice, J; the latter are also offset by half a lattice spacing in the 
direction out of the paper (the /i = r direction). The all zero for 

H = T,x, while the only non-zero values of ajy are shown above. Notice that 
the a° flux obeys (ITT|1 . 

© by the Poisson summation formula: 




^ef,^xaJ^l^uAjx ,(10) 



□ 

where aj^ (like qj^) is an integer-valued vector field on the links of the dual 
lattice (here, and below, we drop overall normalization factors in front of 
the partition function). Next, we write the Berry phase in a form more 
amenable to duality transformations. Choose a 'background' aj^ = flux 
which satisfies 

where j is the direct lattice site in the center of the plaquette defined by 
the curl on the left-hand-side. Any integer-valued solution of (fTT|) is an 
acceptable choice for Oj^, and a convenient choice is shown in Fig^ Using 
f|TT|l to rewrite the Berry phase in ©, applying (fTIH) . and shifting a^^ by the 



6 



• • • • • 

1/4 1/4 

• • • • • 

3/4 1/2 3/4 1/2 

• • • • • 
1/4 1/4 

• • • • • 

(a) 

Figure 2: Specification of tlie non-zero values of tlie fixed fields (a) Xj and 
(b) introduced in ()14|) . The notational conventions are as in Fig^ Only 
tlie n = T components of are non-zero, and these are shown in (b). 

integer 2Saj^, we obtain a new exact representation of Z in Q: 

(12) 

The integral over the Aj^ can be performed independently on each link, and 
its only consequence is the imposition of the constraint e^y\/S.yaj\ = 0. We 
solve this constraint by writing aj^ as the gradient of a integer- valued 'height' 
hj on the sites of the dual lattice, and so obtain 

Z = J2^W (-^$^(AA-2^ayM . (13) 

{hj} \ j,ti J 

This is the promised 2+1 dimensional interface, or height, model in almost 
its final form. 

The physical properties of (fT^ become clearer by converting the "frus- 
tration" a?^ in (jl3|) into offsets for the allowed height values. This is done by 
decomposing a?^ into curl and divergence free parts and writing it in terms 
of new fixed fields, and as follows: 

= A^Xj + e^,xA,yjx. (14) 
The values of these new fields are shown in Fig El Inserting ()14|) into ()13|). 



+1/8 -1/8 +1/8 -1/8 
XXX 

-1/8 +1/8 -1/8 +1/8 
XXX 

+1/8 -1/8 +1/8 -1/8 
(b) 
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we can now write the height model in its simplest form jj] 



Zh = ^ exp 



(15) 



where 



H: 




(16) 



is the new height variable we shall work with. Notice that the have 
dropped out, while the Xj act only as fractional offsets (for S not an even 
integer) to the integer heights. From (fT^ we see that for half-odd-integer 
S the height is restricted to be an integer on one of the four sublattices, 
an integer plus 1/4 on the second, an integer plus 1/2 on the third, and 
an integer plus 3/4 on the fourth; the fractional parts of these heights are 
as shown in Fig the steps between neighboring heights are always an 
integer plus 1/4, or an integer plus 3/4. For S an odd integer, the heights 
are integers on one square sublattice, and half-odd-integers on the second 
sublattice. Finally for even integer S the offset has no effect and the height 
is an integer on all sites. We discuss these classes of S values in turn in the 
following subsections. 

3.1 S even integer 

In this case the offsets 2SXj are all integers, and (jisp is just an ordinary three 
dimensional height model which has been much studied in the literature. 
Unlike the two-dimensional case, three-dimensional height models generically 
have no roughening transition, and the interface is always smooth. With all 
heights integers, the smooth phase breaks no lattice symmetries. So square 
lattice antiferromagnets with 5" even integer can have a paramagnetic ground 
state with a spin gap and no broken symmetries. This is in accord with the 
exact ground state for a 5 = 2 antiferromagnet on the square lattice found 
by Affleck et ai, the AKLT state 0. 



Now the heights of the interface model can take four possible values, which 
are integers plus the offsets on the four square sublattices shown in Fig|2ti'- 
As in Section 13. H the interface is always smooth i. e. any state of (fT3j) has 



3.2 S half-odd-integer 
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Figure 3: Mapping between the quantum dimer model and the interface 
model Z in (fT3j) . Each dimer on the direct lattice is associated with a step 
in height of ±3/4 on the link of the dual lattice that crosses it. All other 
height steps are ±1/4. Each dimer represents a singlet valence between the 
sites, as in Fig|21 

a fixed average interface height ^^{H^-, and any well-defined value for this 
average height breaks the uniform shift symmetry of the height model under 
which ifj — > ifj ± 1. After accounting for the height offsets, we will see 
below that any smooth interface must also break a lattice symmetry with 
the development of hond order, this allows a number of distinct spin gap 
ground states of the lattice antiferromagnet. 

It is useful, first, to obtain a simple physical interpretation of the interface 
model in the language of the 5 = 1/2 antiferromagnet [7j. From Fig|2ti' it is 
clear that nearest neighbor heights can differ either by 1/4 or 3/4 (modulo 
integers). To minimize the action in (jl5|) . we should choose the interface 
with the largest possible number of steps of ±1/4. However, the interface is 
frustrated, and it is not possible to make all steps ±1/4 and at least a quarter 
of the steps must be ±3/4. Indeed, there is a precise one-to-one mapping 
between interfaces with the minimal number of ±3/4 steps (we regard inter- 
faces differing by a uniform integer shift in all heights as equivalent) and the 
dimer coverings of the square lattice: the proof of this claim is illustrated in 
Fig 01 We identify each dimer with a singlet valence bond between the spins 
(the ellipses in Fig 12)), and so each interface corresponds to a quantum state 
with each spin locked in the a singlet valence bond with a particular near- 
est neighbor. Fluctuations of the interface in imaginary time between such 
configurations correspond to quantum tunneling events between such dimer 
states, and an effective Hamiltonian for this is provided by the quantum 
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(a) (b) 

Figure 4: Sketch of the two simplest possible states with bond order for 
5 = 1/2 on the square lattice: (a) the columnar spin-Peierls states, and 
(b) plaquette state. The different values of the (Sj • Sj) on the links are 
encoded by the different line styles. Both states are 4-fold degenerate; an 
8-fold degenerate state, with superposition of the above orders, also appears 
as a possible ground state of the generalized interface model. 

dimer model [Hj. 

The nature of the possible smooth phases of the interface model are easy 
to determine from the above picture and by standard techniques from sta- 
tistical theory lUEI. Interfaces with average height {Hj) = 1/8,3/8,5/8,7/8 
(modulo integers) correspond to the four-fold degenerate bond-ordered states 
in FigE^, while those with (Hj) = 0,1/4,1/2,3/4 (modulo integers) corre- 
spond to the four-fold degenerate plaquette bond-ordered states in FigHjD. All 
other values of {Hj) are associated with eight-fold degenerate bond-ordered 
states with a superposition of the orders in Fig|^ and b. 

Support for the class of bond-ordered states described above has appeared 
in a number of numerical studies of S = 1/2 antiferromagnets in d = 2 
which have succeeded in moving from the small g Neel phase to the large 
g paramagnet. These include studies on the honeycomb lattice P|, on the 
planar pyrochlore lattice ^Hj, on square lattice models with ring-exchange 
and easy-plane spin symmetry [T2], and square lattice models with SU(iV) 
symmetry |14j . 

3.3 S odd integer 

This case is similar to that S half-odd-integer, and we will not consider it in 
detail. The Berry phases again induce bond order in the spin gap state, but 
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this order need only lead to a two-fold degeneracy. 



4 Conclusions 

The primary topic discussed in this paper has been the effective field the- 
ory of paramagnetic Mott insulators with coUinear spin correlations. This 
field theory is the compact U(l) gauge theory in (jHI), and applies in all spa- 
tial dimensions. We also reviewed duality mappings of (jH)) which are spe- 
cial to d = 2 spatial dimensions, and mapped the theory onto the interface 
model ()15p. Finally, we reiterate that paramagnetic Mott insulators with 
non-coUinear spin correlations are described by a Z2 gauge theory which has 
not been presented here. 
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